Flexive and Propulsive Dynamics of Elastica at Low Reynolds Numbers 
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A stiff one-armed swimmer in glycerine goes nowhere, but if its arm is elastic, exerting a restorative 
torque proportional to local curvature, the swimmer can go on its way. Considering this happy 
consequence and the principles of elasticity, we study a hyperdiffusion equation for the shape of 
the elastica in viscous flow, find solutions for impulsive or oscillatory forcing, and elucidate relevant 
aspects of propulsion. These results have application in a variety of physical and biological contexts, 
from dynamic biopolymer bending experiments to instabilities of elastic filaments. 

PACS numbers: 03.40.Dz, 47.15.Gf, 87.45.-k, 
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In Stokes flow, the Aristotelian fluid regime inhabited 
by the very small or the very slow, inertia is irrelevant. 
This fact underlies the inability of a variety of swimming 
motions, perfectly successful on human scales, to gener- 
ate net motion on microscopic scales [0. An oft-quoted 
example is the lack of propulsion for a swimmer with only 
one degree of mechanical freedom, e.g. the paradigmatic 
scallop of Purcell's 1977 lecture which introduced many 
to the principles of Stokes flow Q. Colloquially known 
as "the scallop theorem," this observation derives from 
the more general statement that motions invariant under 
t — ► —t can produce no net effect pj; movies of Stokes 
flow must appear equally sensible when reversed H. 

Purcell observed two ways to elude the scallop theo- 
rem: rotate a chiral arm, or wave an elastic arm. While 
the former dynamic is well-studied (most notably, in the 
context of E. coli [Q ) , the latter is largely uninvestigated 
f|, despite its relation to experiments from motility as- 
says to dynamic studies of biopolymer bending moduli. 
To elucidate this dynamic, we here study the motion 
of a one-armed swimmer with an elastic prosthesis, or 
equivalently a driven elastic filament. Building on exper- 
iments showing the hyperdiffusivity of small-amplitude 
planar deformations p], we quantify how an elastic fil- 
ament eludes the scallop theorem, suggest experiments 
to test these results, and show how this analysis allows 
for the measurment of bending moduli. Remarkably, the 
required mathematics j?]] is central to an intrinsic de- 
scription of overdamped elastica in three dimensions. 

Force- velocity proportionalities in Stokes flow are gen- 
erally not simple; notable exceptions are for highly sym- 
metric objects such as the sphere (F — 6nfiav, with a 
the radius) , and those for which length L greatly exceeds 
width d, where slender-body hydrodynamics [^|j9| applies. 
To lowest order in l/\n(L/d), force and velocity obey 
a local, anisotropic proportionality. For velocity v nor- 
mal to the long axis, the force per unit length / = Qv, 
with ( = 47r/j,/(ln(L/d) + c), where c is an 0(1) con- 
stant determined by the shape of the object. For an 
clastic filament with bending modulus A, —/is the func- 



.2 

XX ' 



tional derivative of the bending energy A/2 J Q dx y 
written here for small planar deformations y{x); thus, 
/ = — Ay xxxx . At free ends the functional derivative 
implies boundary conditions of torquelessness and force- 
lessness: y xx = y xxx = 0. For small deformations v = y%, 
and with the hyperdiffusion constant v — A/C,, we have 



yt = -vy. 



XXXX ? 



(1) 



perhaps the simplest model for the balance between vis- 
cous drag and bending elasticity. 

In 1851, Stokes suggested two problems in fluid me- 
chanics, here termed SI and SII (Fig. to illustrate 
viscous diffusion of velocity |l0| ] : SI - impulsively move a 
wall bounding a fluid; SII - oscillate the wall at frequency 
lj. These motivate two problems for the elastohydrody- 
namic equation of motion (|l|) : EHDI - impulsively move 
one end of a filament; EHDII - oscillate the end. In SI 
and SII, the Navier-Stokes equation reduces to a diffusion 
equation ut — vu xx , with kinematic viscosity v = jij 'p. 
For a semi-infinite domain, the post-transient solution of 
SII consists of decaying, right-moving waves, u (x, i) = 
C/exp(— r]/^/2) cos (t?/v2 — ut), with r\ — xjt v , and vis- 
cous penetration length l v — (v/oj) 1 / 2 . In EHDII the 
analogous elastohydrodynamic penetration length is 

(u/uj) . 



y t = -vy. 



(2) 
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FIG. 1. Geometry of Stokes problems I and II (left) and of 
elastohydrodynamic problems I and II (right). 
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Imposing the left filament end position yo cos tot and 
torquelessness for the left end [[yj y xx (0, i) = 0, we find 
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e dr > cos ( Srj + wtj + c §r > cos (Crj - wt 



(3) 

where C — cos(7r/8), S = sin (71-/8), and now r\ — xji^- 
Unlike SII, EHDII supports left- and right-moving waves 
(with different velocities and decay lengths), despite the 
lack of a reflecting right-end boundary. 

For finite filaments we define a rescaled length £ = 
L/£p(uj) and coordinate a = x/L. When £ < 1 the 
filament behaves as a rigid rod, while it undulates appre- 
ciably for £ ^ f . In this way £c resembles the persistence 
length L p . The exact solution of EHDII for finite L |7|] 
has an expansion in powers of £ whose first terms are 



JJ_ 

yo 



3 \ £4 

—a ) cos(uit) H 1 

2 j V ' 1680 



16a - 70a d 



-70a 4 



21a 5 ) sin(wt) + 0(£ 8 ) . (4) 

At order £°, the filament is a straight rod that interest- 
ingly pivots about a point two-thirds of its length. Flex- 
ive corrections at 0(£ 4 ) break time-reversal invariance. 
Solutions of increasing £ are shown in Fig. |^, whose 
inset shows the results of an experiment on actin || in 
which observed shapes in a range of frequencies were fit 
for ly to the exact expressions 0, verifying the scaling 
l„ ~ lu^ 1 / 4 as well as providing a novel dynamic tech- 
nique for measuring the bending modulus A. 

The propulsive force F x imparted to the fluid by the 
filament (or vice versa) may be computed by integrating 
the projected elastic force density f along the filament. 
Interestingly, with the boundary conditions of EHDII and 
the geometrically exact form of f this result is expressible 
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FIG. 2. Solutions to (|l|) for oscillatory driving at various 
rescaled lengths C = L/l„. Inset: experimental verification 
M of the frequency dependence of £h(uj) yieldinga modulus 
/ksT — 7.4/im, well within the range (0.5/im |12] to 17/im 
of measurements employing statistical techniques. 



in terms of the curvature k and tangent angle 9 at the 
forcing point: F x — Ak s sui9(s = 0) ~ Ay x y xxx (x = 0). 
The time average of this quantity over one period gives 



F x = 



■Y(C) 



(5) 



For £ < 1, Y ~ (ll/3360)£ 4 , so F x - y%( 2 Lu 2 L 4 /A: a 
short (or infinitely stiff) pivoting filament produces no 
net force (by the scallop theorem). Flexibility leads to a 
net leftward propulsion, as the right-moving waves domi- 
nate the left-movers. An unexpected, fascinating feature 
(Fig. |3|) is the maximum in Y indicating an optimal value 
of the length L* ~ 4.07 (A/ (u) 1 / 4 . 

In a familiar way, this force is associated with the 
trajectory of the filament shape in a low-dimensional 
projection of configuration space. The relation F(t) = 
Ay x (0,t)y xxx (0,t) and the equation of motion imply 



F x = 2- 
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where A is the area under the curve y{x) and 9q ~ y x 
is the tangent angle at the left end. Thus, the net force 
during the cyclic motion is the area enclosed in the asso- 
ciated 'Carnot diagram' in A — 9q space; it results from 
pushing aside some volume of fluid (area, in two dimen- 
sions), projected via #0 in the direction of propulsion. For 
EHDII, the trajectory is an ellipse; it thins to a straight 
line for the time-reversible pivoting of a rod, encloses no 
area, and thus produces no force. Observe the intuitive 
result that net propulsion in the transverse direction, pro- 
portional to § dA = 0, vanishes identically. 

We estimate the efficiency £ of this motion Q by com- 
paring the power Pu = F x v x ~ F 2 /LC,\\ for longitudinal 



propulsion to the power P± 
pated in transverse motions, to obtain 
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where Z(C) is the scaling function shown in the inset 
to Fig. Q. Filaments short relative to £v(u>) flex little 
and produce little propulsion, while long ones have excess 
drag from the nearly straight regions far away from the 
point of forcing, yielding a sharp maximum at £ ~ 4.0. 

These observations suggest experiments in the spirit of 
those by Taylor on swimming with a helical flagel- 
lum. Exploiting the results of EHDII, perhaps carried 
out on microfilaments with laser || or magnetic tweez- 
ers, or on macroscopic objects, one might measure the 
propulsive force through the transverse displacement at 
the forcing point, test for the predicted maximum as a 
function of frequency, investigate the role of nonlinear- 
ities, and study interactions between flexing filaments. 
Analogous experiments incorporating twist (perhaps via 
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magnetic optically-trapped beads, as in |Tf|) could elu- 
cidate instabilities exhibited by helical motion of flexible 
filaments p3] and associated propulsion. 
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FIG. 3. Scaling functions Y, for propulsive force, and Z, 
for efficiency, versus rescaled length £ — L /£„(&). 

A swimmer's arm may also move by discrete, noncyclic 
strokes. The associated model dynamic (EHDI) is the re- 
laxation of a filament to an equilibrium shape given some 
boundary condition at the point of attachment. Since 
(Q) is linear, we may subtract from its general solution 
a particular solution consistent with nonzero boundary 
conditions or external driving to obtain a homogeneous 
problem with 0-valued boundary conditions. This mo- 
tivates the construction of a self-adjoint operator from 
Ti = d x , whose well-known eigenfunctions are 

yV q (x) — a\ sin (qx) + a 2 cos (qx) 

+a% sinh (qx) + 04 cosh (qx) , (8) 

where HW q = q 4 W q . The ten distinct choices of bound- 
ary conditions for which Tt is self adjoint deter- 
mine the distinct coefficients {a^} of W q and the allowed 
wavenumbers q. These q are roots of transendental solv- 
ability conditions; certain boundary conditions are satis- 
fied by Fourier modes (those for which 03 = 04 = 0). 

The completeness of the set W q yields the homoge- 
neous solution to (|l]), 

y(x,t) = J2w q (x)e-^ 4t / dxW q (x)y(x,0) , (9) 

relevant to a second class of techniques for measuring 
bending moduli. One recent example fT^ | is the relax- 
ation of the free end of an initially bent microtubule 
whose opposite end is clamped to a cover slip via the 
axoneme from which it is nucleated. Using typical mate- 
rial parameters for microtubules, (L ~ /im, £ ~ cP, L p ~ 



mm), we observe from (g) that beyond the inital few mil- 
liseconds the filament straightens with a time constant 
Ti = \fvq\\ here q\ ~ 1.875/1/ the smallest root of the 
clamped- free solvability condition cos(qi) cosh(qi) = — 1 
p8[ . Measurement of t\ gives A. 

Finally, we apply these results to an important subject 
in pattern formation: intrinsic formulations of filament 
motion p9[ and their extention to nonplanar dynamics. 
The dynamical evolution of the Frenet-Serret (FS) cur- 
vature k and torsion r are singular at inflection points 
(where n = and r is undefined) , a problem removed by 
instead using the parameterization p0|j2l[ | 

tp(s, t) = k(s, i)e^ s '*> , w = (n + ib)e*+ , (10) 

where 0(s, t) = J s c1s't(s' , t), and n and b are the normal 
and binormal vectors. Across an inflection point at, say, 
s = 0, the normal and binormal vectors n and b flip sign, 
the torsion has a singular piece, r = tt5(s), so there is 
a discontinuity <fi(0 + ) — (f>(0~) = tt. Yet, this leaves w 
and tp continuous. The curve r(s) is constructed from tp 
via new FS equations: w s = — tpt; t s = Re("0*w). Some 
elementary shapes have very simple ^-representations: 
the straight line ■0 = 0, the circle ip = a, the helix ip = 
ae lps (a,p constant). 

The intrinsic inextensible evolution for ip po| , 

rlH = (d ss + iV-f )T + V-dm J ds'^T* + W^ s , (11) 

depends upon the velocity components Yt — Uh + Vh + 
Wt through the complex velocity T = (U + iV)e 1 ^. It 
can be shown p2| from the governing equations for vis- 
cous motion of filaments, that T = —i>(tp ss + (1/2)1^ \ 2 ip). 
Thus, the ip dynamics is hyperdiffusive as in (|l|) [ p0| | , 

ipt ^ -i>4> SSS s H • (12) 

For an elastic energy £ = (A/2) J dsn 2 the boundary con- 
ditions at free ends are k = k s = kt = 0, or simply 
ip = ip s = 0. Since the self-adjoint operator is again 
H = dg, the W q (s) 1 now functions of arclength, are the 
relevant eigenfunctions. Thus, the i[> formulation, the 
natural singularity-free mathematical representation of 
the curve, also compactly encodes the physical boundary 
conditions and dynamics. Since the W q form a complete 
set, the evolving shape of any filament is expressible as 

i)( S ,t) =Y,C q (t)Wg(s) , 
9 

the time evolution being a nonlinear dynamical sys- 
tem in c q (t) = (yV q \?p(s, t)); in the linearized regime, 
c q (t) — c q (0) exp(— vq A t). The modes W q are like the 
Rouse modes of a flexible polymer, and the "clamped" 
eigenfunctions of yt = ~vy xxxx (with y = y x = at 
x = 0,L) are the "free" eigenfunctions of ipt = —vipssss. 
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with cos(qL) cosh(qi) = 1. The first three are shown in 
Fig. ||. When the complex coefficients c q share a com- 
mon phase the curve is planar, otherwise it has torsion. 
The latter case is shown in Fig. | as a coiled three- 
dimensional filament composed of the first two modes, 
reconstructed from the complex FS equations. The sim- 
plification achieved with the ip formulation is clear. 
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FIG. 4. Biharmonic eigenfunctions and a nonplanar elas- 



tica. Top: W 9 with boundary conditions W = d a W 
Bottom: free filament with ip — 15(W gi + iW q2 ). 
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This formalism easily generalizes to twisted filaments 
p2| , as found for instance in bacterial systems that un- 
dergo hierarchical buckling and writhing instabilities |p3[ . 
For a constant twist density Q and twist elastic constant 
C the moment relation M = Anh + Cflt can be shown 
to yield the complex velocity p4fl 



T 



1 



(13) 



so the intrinsic dynamics ( |12| ) is explicitly complex 

ipt ~ -i>(tp sass + ifiP0 sss ) H , 



(14) 



with r = C I A. A helical perturbation ip ~ exp(iks + at) 
to a straight filament leads to a growth rate a = —v{k A + 
Tilfc 3 ). This describes a writhing instability on a length 
scale Lwr = 27r/f2r to a helix whose handedness is set 
by the sign of the twist density f2. A stability analysis 
p2[ about a loop of radius R yields not only the criterion 
TflR = %/3 for the onset of the primary instability p5[ , 
but also the growth rates o~± for the coupled bending and 
writhing modes. Extensions to the fully nonlinear regime 
are then straightforward, leading to a "dynamics of twist 
and writhe" that complements important existing Hamil- 
tonian formulations |26|. 
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